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Abstract 



■ Since the seminal work [19, 20] the iteration stable (STIT) tessellations have 
Ph . attracted considerable interest in stochastic geometry as a natural and flexible, yet 

analytically tractable model for hierarchical spatial cell-splitting and crack-formation 
processes considered in stochastic geometry. We provide in this paper a fundamen- 
tal link between typical characteristics of STIT tessellations and those of suitable 
mixtures of Poisson hyperplane tessellations using martingale techniques and general 
theory of pure jump Markov processes. As applications, new mean values and new 
. distributional results for STIT tessellations are obtained. 
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1 Introduction 

Infinite divisibility or stochastic stability of a random object under a certain operation is 
one of the most fundamental concepts in probability theory, prominent examples include 
the classical theory of infinite divisible and stable distributions with their applications 
around the central limit theorem, max-stable distributions studied in extreme value theory 
or union infinitely divisible random sets studied in the classical theory of random closed 
sets [14]. 

In the present paper we deal with a class of non- stationary iteration infinitely divisible 
random tessellations of the d- dimensional Euclidean space and, more specifically, with sta- 
tionary random tessellations that are stable under the operation of iteration - called STIT 
tessellations for short. Recall that a tessellation (or mosaic) of M> d is a locally finite familiy 
of compact convex random polytopes, which have pairwise no common interior points, 
cover the whole space and that they are a central object studied in stochastic geometry 
and related fields, see [1, 2, 4, 6, 7, 8, 9] to name just a few. 

In the stationary case, the motivation for random STIT tessellations goes back to R. Am- 
bartzumian in the 80thies and the principle of iteration of tessellations can roughly be 
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explained as follows: Take a random primary or frame tessellation and associate with each 
of its cells independently a copy of the primary tessellation, called component tessellation, 
which is also independent of the primary tessellation as well. Perform now in each cell a 
local superposition of the primary tessellation and the associated tessellation, whereupon 
scale the resulting random tessellation by a factor 2 in order to ensure that the mean sur- 
face measure of cell boundaries stays constant. The described operation can be repeatedly 
applied and we obtain in this way a sequence of random tessellations. It can be shown that 
this sequence converges to a random limit tessellation and that this tessellation is stable 
under iterations - a STIT tessellation - in the translation-invariant set-up. In the general 
case, the resulting tessellation is iteration infinitely divisible in any finite volume. 
Starting with [19], STIT tessellations and their theoretical framework were formally in- 
troduced in [20] by Nagel and Weiss. The same research group discovered a first basic 
technique for studying mean values and even some distributions related to the geometry of 
STIT tessellations by writing certain balance equations based on the stochastic stability of 
the tessellation, see [21, 22, 32] for the mean values as well as [15, 16] for distributional re- 
sults. In [17, 18] a new aspect was introduced into the theory, namely a tessellation- valued 
random Markov process on the positive real half-axis with the property that at each time 
the law of the tessellation is stable under iteration. In finite compact convex volumes 
W C M. d with VoLj(W) > 0, this process can be explained as follows: Let us fix a terminal 
time t > and a (in some sense non-degenerate) translation-invariant measure A on the 
space of hyperplanes in R d . In a first step, we assign to W a random lifetime. Upon expiry 
of its lifetime, the primordial cell W dies and splits into two polyhedral sub-cells separated 
by a hyperplane hitting W, which is chosen according to the normalized distribution A. 
The resulting new cells are again assigned independent random lifetimes and the entire 
construction continues recursively until the deterministic time threshold t is reached (see 
Figure 1 for an illustration). The resulting random tessellation is denoted by Y(tA,W). 
In order to ensure the Markov property of the above construction in the continuous-time 
parameter t and in order to have distributional equality of the outcome with the above 
mentioned limit tessellations, we will have to take care of the special choice of the lifetime 
distributions, see the details in Section 3 below. 

This dynamic point of view can be exploited to establish new results on STIT tessellations 
and in fact our theory will be based on it. In addition, the finite volume Markovian con- 
struction provides a link to a class of more general and non- stationary iteration infinitely 
divisible random MNW-tessellations, where the translation-invariant hyperplane measure 
A from above is replaced by an arbitrary non-atomic and locally finite measure on the 
space of hyperplanes. We would like to emphasize that the dynamical representation is 
a special feature of STIT tessellations and their infinitely divisible counterparts and, as 
recently pointed out by the first author [25], that such and similar spatio-temporal random 
processes have remarkable potential for applications in stochastic geometry. 
The purpose of the present paper is to explore the dynamic representation of STIT tes- 
sellations further and to introduce a new technique, which has has the advantage to deal 
with properties of the model that were out of reach so far. The crucial observation is 
that the dynamic representation of the tessellations under consideration fulfils a Markov 
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property in the continuous-time parameter, which allows to construct certain martingales 
related to the tessellations (see Section 4), eventually leading to fundamental comparison 
results of STIT tessellations with certain mixtures of Poisson hyperplane tessellations in 
Section 5. As an application of these results we calculate new mean values and also some 
new distributions related to geometric objects determined by the tessellation. To keep the 
paper self-contained, we recall in Section 2 the construction of STIT tessellations as limit 
of repeated iterations of tessellations and formally introduce in Section 3 their Markovian 
dynamic representation. Moreover, we also introduce there the class of iteration infinitely 
divisible random MNW-tessellations and summarize some of their most important proper- 
ties needed in this paper. 

The current work is based on an extended version available online [26]. It also forms the 
basis of our subsequent papers [27, 28, 29, 30]. 

2 STIT tessellations as limits 

A tessellation of IR d is a locally finite partition of the space into compact convex polytopes, 
the cells of the tessellation. One can regard a tessellation either as a collection of its cells 
or as the closed set formed by the union of their boundaries. We will mostly follow the 
second mentioned point of view and denote by Cells(Y) the set of cells of the tessellation Y 
(note that by the Jordan-Schonflies theorem the correspondence between Y and Cells (Y) 
is one-to-one). Thus, a random tessellation can be regarded as a special random closed set 
in the classical sense of stochastic geometry, see [14, 23, 24]. In particular, this imposes 
the usual Fell topology and the corresponding Borel measurable structure on the family 
of tessellations, see ibidem. A random tessellation Y (regarded as a random closed set in 
M. d ) is called stationary if its distribution does not change upon actions of translations. 
Analogously a random tessellation is called isotropic if its distribution is invariant under 
the action of SO(d). 

Whenever we two random tessellations Y 1 and Y 2 of M. d are given we can define their 
iteration/nesting. To this end, we associate to each cell c G Cells(Yi) an independent copy 
Y 2 (c) of Y 2 and we assume furthermore the family {Y 2 (c) : c £ Cells(Yi)} to be independent 
of Y\. Then we define the iteration of Y\ with Y 2 by 

Yi mY 2 :=Y 1 m{Y 2 (c) : cG Cells(Yi)} :=Y 1 U \J (Y 2 (c)nc), 

ceCells(Yi) 

i.e. we take the local superposition of Y\ and the family {Y 2 (c) : c G Cells(Yi)} inside the 
cells of Y\. It was shown in [18] that Y\ EBY2 is a stationary random tessellation as soon Y\ 
and Y 2 have this property. Note also that the so-defined nesting operation is associative in 
distribution. A stationary random tessellation Y is called stable under iteration, or STIT 
for short, iff 

m(Ya...SY) =Y, to = 2, 3,..., (1) 

v v y 

m 
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where = stands for equality in distribution (note that rescaling by factor m ensures that 
the mean surface area of cell boundaries per unit volume remains constant). In fact, using 
the uniqueness results, see Theorem 3 and Corollary 2 in [20], it is easy to see that it is 
enough to take one fixed m > 1 in (1). 

To proceed, let us be given a constant < t < oo and a probability measure 1Z on the unit 
sphere Sj-i, usually identified with the induced distribution of orthogonal hyperplanes on 
the space of (d— l)-dimensional linear hyperplanes in M. d , also denoted by 7Z in the sequel 
for notational simplicity. Define the measure A on the space H of afline hyperplanes in M. d 
as the product measure 

A:=£ + ®TZ (2) 

of £ + standing for the Lebesgue measure on the positive real half-axis (0, oo), and of 7Z, 
where a pair (r,u) G (0, oo) x Sd-i is identified with the hyperplane {x G M. d , (x,u) = r}. 
Throughout this paper we always require that the support of 1Z spans the whole space, i.e. 
that span(supp(7?.)) = M. d . Assume now that we are given a stationary random tessellation 
Y with surface intensity t (i.e. the mean surface area of cell boundaries per unit volume 
equals t) and directional distribution 1Z (i.e. the distribution of the normal direction of 
the facet containing the typical point in the standard Palm sense is given by 1Z) and define 
the sequence (2 n (y)) n > 1 by 

Ti 

XAY) := 2(Y ffl Y), l n (Y) := l n -i(Y) ffl nY = n (Y ffl . . . ffl Y), n>2. 

n — 1 » v ' 

n 

It was shown in [20, Thm. 3] that X n {Y) converges in law, as n — > oo, to a stationary 
random limit tessellation Y(tA) uniquely determinedby tA. This tessellation is easily shown 
to be stable under iterations, whence a STIT tessellation with parameters t and A and we 
refer to Figure 1 for an illustration. 

3 The MNW-construction 

As already emphasized in the Introduction, it is a crucial feature of the random tessellations 
Y(tA) introduced in the previous section that they admit a natural and intuitive spatio- 
temporal construction. For a restriction Y(tA, W) of Y(tA) to a compact convex window 
W C M, d with interior points, this construction can be described as follows (the reader 
is referred to [20] for full details): Assign to the window W an exponentially distributed 
random lifetime with parameter A([W]) where [W] := {H G "H, H D W ^ 0} stands for 
the family of all hyperplanes hitting W . Upon expiry of its lifetime, the cell W dies and 
splits into two sub-cells W + and W~ separated by a hyperplane in [W] chosen according to 
the law A(-)/A([W]). The resulting new cells W + and W~ are again assigned independent 
exponential lifetimes with respective parameters A([IU + ]) and A([VF - ]) (whence smaller 
cells live stochastically longer) and the entire construction continues recursively until the 
deterministic time threshold t is reached. The cell-separating (d— l)-dimensional facets (the 
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Figure 1: Realizations of a planar and a spatial stationary and isotropic STIT tessellation 
(kindly provided by Joachim Ohser and Claudia Redenbach) 



word facet stands for a (d — l)-dimensional face here and throughout) arising in subsequent 
splits are usually referred to as (d — 1) -dimensional maximal polytopes (or I-segments for 
d = 2 as assuming shapes similar to the letter I). The described process of recursive 
cell divisions is called the Mecke-Nagel- Weiss- or MNW- construction in the sequel and the 
resulting random tessellation created inside W is denoted by Y(tA, W) as mentioned above, 
whereas the collection of all (d— l)-dimensional maximal polytopes or I-segments is denoted 
by MaxPolytopes d _ 1 (F(tA, W)). Moreover, we write MaxPolytopes fc (F(tA, W)) for the 
collection of k- dimensional maximal polytopes of Y(tA,W), where by a A> dimensional 
maximal polytope we mean the maximal union of connected and A;-coplanar k- dimensional 
faces of cells. In fact, the A;- dimensional maximal polytopes of Y(tA, W) are nothing than 
the fc-faces of (d — l)-dimensional maximal polytopes. 

It was shown in [20] that the law of Y(tA, W) is consistent in that Y(tA, W) n V = Y(tA, V) 
for convex V C W and thus Y(tA, W) can be extended to random tessellation Y(tA) on 
the whole space, which is then proved (see [20]) to coincide with the limit tessellation 
Y(tA) considered in the previous section, as the notation suggests. Again, the sets of 
all (d — l)-dimensional and fc-dimensional maximal polytopes of Y(tA) are denoted by 
MaxPolytopes d _ 1 (F(tA)) and MaxPolytopes fc (F(tA)) (0 < k < d - 2), respectively. The 
stationary random tessellation Y(tA) is additionally isotropic if and only if 1Z is the uniform 
distribution on Sd~i in the factorization (2). 

A simple yet crucial observation is that even though only translation-invariant measures A 
of the form (2) show up in the limiting STIT tessellations, the MNW-construction can be 
carried out with arbitrary non-atomic and locally finite driving measure A (i.e. A([W / ]) < oo 
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for W C M. d bounded) on H also leading to a consistent family Y(th, W) and eventually, by 
extension, yielding Y(tA). Many of our theorems will be stated in this general context. It 
should be emphasized though that such tessellations are no longer iteration stable (STIT). 
However, they have the general property of being iteration infinitely divisible, as they can 
be readily checked to arise as m-fold iterations of Y(t/m) for each m > 2 in all finite 
volumes. Formally, this means that 

Y((t/m)A,W) mm = Y(tA,W) 

for all compact convex sets W C M d , which follows directly by the MNW-construction as 
yielding 

Y(sA, W) ffl Y(uA, W) = Y((s + u)A, W). 

It is more than natural to expect that also Y(tA) = Y((t/m)A) mm for all m in the whole 
of M. d , which should be provable by adopting the theory developed in [18], thus even better 
justifying the term iteration infinitely divisible in our context, yet this falls beyond the 
scope of the present work and is not needed in our arguments. 

It is worth pointing out that it is currently an open problem whether any iteration in- 
finitely divisible random tessellation in R d can be constructed with the MNW-algorithm. 
This means that in this paper infinite divisibility with respect to iteration only refers to 
such tessellations which can be obtained by the MNW-construction. For this reason, the 
tessellations we are dealing with in the general set-up are called iteration infinitely divisible 
MNW-tessellations. 

4 Associated martingales 

The finite volume continuous-time incremental MNW-construction of iteration infinitely 
divisible random tessellations or more specially stationary STIT tessellations, as discussed 
in Section 3 above, clearly enjoys the Markov property in the time parameter, whence 
natural martingales arise, which will be of crucial importance for our further considerations. 
To discuss these processes we need some additional terminology and notational conventions. 
Let us fix a compact convex window W and a general diffuse (non-atomic) and locally finite 
measure A on H, the space of hyperplanes in M. d . Next, for a tessellation Y of M. d and any 
H G [W], we denote by Cells(F D H) the set of cells the hyperplane H is tessellated with 
as a result of the intersection with Y. In terms of the MNW-construction, the members 
of Cells(y fl H), which are (d — l)-dimensional polytopes, are the potentially new facets 
to be added to the tessellation when the construction is continued. Whenever a new facet 
/ G Cells(y H H) to be added to Y should a facet birth (cell split) occur on H in Y, we 
denote by Cell(/, H\Y) the cell splitting by the birth of / and we write Cell + (/, H\Y) and 
Cell~(/, H\Y) for the two sub-cells into which Cell(/, H\Y) divides, lying on the positive 
and negative sides of H, respectively. 

With this notation, it is readily seen that for a fixed measure A, (Y(tA, W)) t >o is a pure 
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jump Markov process with values in the space of tessellations in W and with the property 
that 



dF (Y((t + dt)A) = YU {f}\Y(tA) =Y) = l[f e Cells(Y n H)]A(dH)dt, (3) 
with H e[W], and 



where | Cells (Y n H)\ is the cardinality of the set Cells (Y n H). Indeed, this is because, 
conditionally on Y(tA) = Y, during the period (t,t + dt] of the MNW-construction we have, 
as discussed in Section 3: 

• For each cell c G Cells(Y) the probability that it undergoes a split is A([c])dt, 
moreover the probability that two or more cells split is o(dt), whence the prob- 



• Should a cell c G Cells (Y) split, the splitting hyperplane H is chosen according to 
the law 1[H G [c]](-)/A([c]), whence the probability of observing a split of c induced 
by H during the time period (t, t + dt] is just 1[H G [c]](dH)dt. Thus, we see that 
on the event Y(tA) = Y, having Y((t + dt)A) = Y U {/} with / C H, H G [c] and 
c G Cells (Y) is equivalent to having c split by H during (t, t + dt]. As noted above, 
the latter happens with probability A(dH)dt whence (3) follows. 

It should be emphasized at this point that the differential notation in (3) and (4) employing 
the symbols cflP and dt is widely accepted in the theory of pure jump continuous-time 
Markov processes and makes perfect formal sense as a commonly recognized abbreviation 
for the usual description in terms of waiting times, with dF(Y((t + dt)A) = Y'\Y(tA) = 
Y) = a(Y, Y')dt understood as : 'while in state Y, wait an exponential time with parameter 
a(Y, Y') and then jump to Y' unless some other jump has occurred prior to that '. We refer 
the reader to Chapter 15 in [5] and especially to Section 15.6 there, where this construction 
is formalized with full mathematical rigour in terms of waiting times as noted above. As 
readily verified, specializing the generic waiting-time construction to the case of (3) and 
(4) yields precisely the standard MNW construction from Section 3. 

Using (3) and (4) we conclude by general theory of Markov processes and their infinitesimal 
generators, see Chapter 1 in [13] or Chapter 15 and especially Sections 15.4 (Def. 15.21) and 
15.6 in [5] specialized for the pure jump case, the generator for (Y(tA, W)) t >o is L := h\. w 
with 






(5) 



7 



for all F bounded and measurable on space of tessellations of W. Consequently, again by 
standard theory as given in Lemma 5.1 Appendix 1 Sec. 5 in [11], see also Section 1.5 
in [13], or alternatively by a direct check straightforward in the present set-up, we readily 
conclude 

Proposition 1 For F bounded and measurable, the stochastic process 

F(Y(tA,W))- [ LF(Y(sA,W))ds (6) 
Jo 

is a martingale with respect to the filtration $t t generated by (Y(sA, W))o< s < t . 
To proceed towards the crucial Proposition 2, consider F of the form 

My) ■= E 

/eMaxPolytopes d _ 1 (Y) 

where, recall, MaxPolytopes (i _ 1 (F) are the (d— l)-dimensional maximal polytopes of Y (the 
I-segments in the two-dimensional case), whereas 0(-) is a generic bounded and measurable 
functional on (d — l)-dimensional facets in W, that is to say a bounded and measurable 
function on the space of closed (d — l)-dimensional polytopes in W, possibly chopped off 
by the boundary of W, with the standard measurable structure inherited from space of 
closed sets in W. Whereas the so-defined F is not bounded and thus (6) cannot be applied 
directly, we can apply it for its truncation F N := (F A N) V — N, N G N, which is 
bounded, and let iV ->■ oo to conclude that F(Y(tA, W)) - f* LF(F(sA, W))ds is a local 
S r martingale, see Definition 5.15 in [12] and take T N = M{t > : \F(Y(tA, W))\ > N} 
there. Now, apply the proof of Lemma 1 in [20], where the number of cells in Y(tA, W), 
and hence for all Y(sA,W), s < t, is bounded by a Furry- Yule-type linear birth process 
whose cardinality at any given finite time admits moments of all orders, to conclude that 
(F(Y(tA, W))- f*LF(Y(sA, W))ds) t < a is of class DL for all a > in the sense of Definition 
4.8 in [12]. Using now Problem 5.19 (i) in [12] we finally conclude that F(Y(tA,W)) - 
J 'LF(F(sA, W))ds is a martingale. Thus, applying (5) and (6) for F = we have 

Proposition 2 The stochastic process 

^(Y(tA,W))- f j £ cf>(f)A(dH)ds 

JO J[W] f eC en S (Y{sA,W)nH) 

is a martingale with respect to 

5 Relationships for intensity measures 

In this section we establish a two fundamental first-order properties of Y(tA, W) for general 
locally finite non-atomic measures A, essentially obtained by comparison with suitable 
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mixtures of Poisson hyperplane tessellations. The key to our results is Proposition 2 from 
Section 4. To exploit it, consider the random measures 



M Y(tA,w). = S c , M Y ^ tA ^ := EM Y{tA ' W) (7) 

ceCells(Y(tA,W)) 

with 8 C standing for the unit mass Dirac measure at c. In full analogy, define - A/f PHT ( tA ' H/ ) 
and M PKT ( tA ' w \ where PHT(tA, W) is the Poisson hyperplane tessellation with intensity 
measure tA, restricted to W (see [23] for background material on Poisson hyperplane 
tessellations). Further, put 

jf W) := ^ ¥ l (tA ' W) ■= ^I itA ' W \ k = l,...,d-l, 

f e MaxPolytopcs k (Y (t A, W)) 

where, recall, MaxPolytopes fc (F) is the collection of /c-dimensional maximal polytopes of 
Y. Likewise, define 

j^m,W) ;= £ F PHT( t A, W ) ;= Ejr PHT( t A )W)) k = li ^ mfd _ lf 

/GFaccs fc (PHT(tA,W / )) 

where Facesfc(PHT(tA, W)) is the collection of all fc-face of the Poisson hyperplane tessel- 
lation PHT(tA, W). Our first claim is 

Theorem 1 We have 

M Y(tA,W) = M PHT(tA,W)_ 

It is worth pointing out that in the translation-invariant set-up, we obtain as a corollary 
that the distribution Q of the typical cell of the STIT tessellation Y(tA) coincides with the 
typical cell distribution Q PHT (* A ) of a stationary Poisson hyperplane tessellation PHT(tA) 
with intensity measure tA, previously shown in [19] by completely different arguments. 
Recall, that the typical cell of a tessellation is a randomly selected cell out of a large 
observation window, where each cell has the same chance of being selected. Indeed, we 
have by the very definition of the distributions Q and Q PHT (* A ) ( c f. [23, Eq. (4.8, 4.9)]) 
and by Campell's theorem [23, Thm. 3.1.2] for any measurable subset A of the space of 
rf-dimensional polytopes 

Acen s(m) )Q(A) = li^ -^-L^ E / l[c — m(c) £ A]M Y ^ tA,rW \dc) 

= lim ^7T7T^ / !t c - m ( c ) G A]M Y ^ rW \dc) 
r->oo r d Vo\ d (W) J 

= lim t- ] — [l[c- m(c) E A]M PUT ^ rW \dc) 
r^oo r d Vo\ d (W) J 1 K ' 

= Acclls(PHT(tA))Q PHT( * A) (-4), 

where Ac e iis(r(tA)) and Ac e iis(PHT(tA)) denote, respectively, the mean number of cells of Y(tA) 
and PHT(tA) per unit volume, W C M. d is measurable with positive volume, and where 
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m(c) stands for a translation-covariant selector for the <i-dimensional polytope c (for ex- 
ample its Steiner point or its center of gravity). Since trivially Ac e iis(r(tA)) — Ac e ii s (PHT(tA)) 
by Theorem 1, we get Q y( * A) = Q PHT (* A ). 

Proof of Theorem 1 Using (5) and (6) with 

F(Y):= ^ 

ceCciis(y) 

for general bounded measurable cell functional 0, with localization argument as the one 
preceding Proposition 2 we conclude that 

1 4>{c')M Y ^ w \dc') -ft Yl mCe\\ + (f,H\Y(tA,W)))+ 

J Jo J[W] f eC ells(Y(sA,W)nH) 

0(Ceir(/, H\Y(tA, W))) - 0(Cell(/, H\Y(tA, W)))]A(dH)ds (8) 

is a ^-martingale. For a polyhedral cell c C W, possibly chopped off by the boundary of 
W, and for H e [c] we write c + (H) and c~(H) to denote the cells into which c gets divided 
by H, lying respectively on the positive and negative side of H. With this notation, (8) 
says that 

/ <P{c')M Y{tA ' w \dc')- f I I [0(c + (iJ)) + 0(c-(/J))-0( C )]A(^)^ y(sA ' w) (^)rfs 

J Jo J J[c] 

is a ^(-martingale. Taking expectations leads to 

,(c')M Y ( tA > w \dc') = III [(f)(c + (H)) + (f)(c-(H))-(f)(c)]A(dH)M Y ^ A ' w \dc)ds 
Jo J J[c] (g) 

for all bounded measurable 0. To proceed, we regard M y ( sA,H/ * ) as an element of the space of 
bounded variation Borel measures on the family of polyhedral sub-cells of W endowed with 
the standard measurable structure inherited from the space of closed sets in W. Consider 
the linear operator Ta on this measure space, given by 

TM =11 [S C+{H ) + S C - [H ) - 5 c ]A(dH)fi(dc). (10) 

J J[c] 

By the definition (10), ||T A (/x)|| TV < (J [w] dA) \\n\\ Ty = A([W]) ||/x|| TV where ||-|| TV is the 
standard total variation norm of a measure, see [3, Def. 3.1.4]. This inequality turns 
into equality when /i — 8w Consequently, Ta is a bounded operator of operator norm 
A([W]) < +oo by the assumed locally finiteness of A. Using the operator T A , the relation 
(9) can be rewritten in form of the operator differential equation 

d M Y(tA,W) = TAM Y(tA,W)^ M Y(0,W) = ^ (n) 
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which, in view of the above properties of 7a, admits by standard theory of linear operators 
(cf. [10, IX.§2, Sec. 2]) the unique solution 

M Y(tA,w) = exp ( tTA } 6wi t > o. (12) 

It is easily seen that exactly the same equations (9), (11) and thus also (12) hold for 
M PHT(tA,wo i n particular, M y (* A - w ) = M PHT (' A < W ) as required. □ 

Having characterized M Y ( tA,w ^ we now turn to ¥ Y ^ tA ' W \ 
Theorem 2 For all k = 1, . . . , d — 1 we have 

t 



¥ Y(tA,W) = {d _ k)2M f l F PHT (s A W(fo 

Jo s 



Proof of Theorem 2 Fix k E {I, . . . ,d — 1}. Let ip be a general bounded measurable 
function of a /c-dimensional maximal polytope, as usual regarded as a closed subset of W, 
and for a (d — l)-dimensional maximal polytope h put 

/€FaceSfc(/i) 

noting that the /c-dimensional maximal polytopes of the tessellation Y(tA, W) are precisely 
the /c-faces of its (d — l)-dimensional maximal polytopes. Using Proposition 2, taking 
expectations and recalling (7) we see that 

E^(Y{tA,W)) = I </)d¥ Y [ tA ' w) = III <f){cnH)A{dH)M Y ^ sA ' w \dc)ds. 

J JO J J[c] 

Applying Theorem 1 we get 

UdW Y[ ^ W) = f (I <f>(cnH)A(dH)M PRT( > sA ' w \dc)ds. (14) 

J Jo J J[c] 

However, by Slivnyak's theory, see e.g. [24, Thm 1.15], we obtain 

J J 0(cn/7)A(rfiJ)M PHT ( A ^)(rfc) = j 0dF™ W) 
and thus, upon taking sA in place of A, more generally, 

J J <p(cnH)A(dH)M PH ^ sA ' w \dc) = -J <f>dF P ^ {sA ' W \ 

whence, with (14), 
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We note now that, by (13), 

J <pd¥ Y d ^ W) = J 4>dF Y k (tA ' W \ 

because each /c-dimensional maximal polytope is a fc-face of precisely one (d— 1) -dimensional 
maximal polytope in Y(tA, W). On the other hand, 



J (f>d¥™ ^ sA ' W) = (d- A;)2 rf - fe - 1 J ^d¥ 



PUT(sA,W) 
k i 



because each £;-face of PHT(sA, W) is a £;-face of (d - k)2 d ~ k ^ facets of PHT(sA, W), see 
Theorems 10.1.2 and 10.3.1 in [23]. Hence, we conclude that 

J ^ ¥ k (tA ' W) = (d- k)2 d ~ k ^ j' lj ^l HT(sA ' W) ds 
for all ip bounded and measurable, which completes the proof of the Theorem. □ 

Some of our arguments in the sequel will require a straightforward formal extension of 
Theorem 2. Namely, we formally mark all (d — l)-dimensional maximal polytopes of the 
tessellation Y(tA, W) by their birth times. This gives rise to the birth-time augmented 
tessellation Y(t,W) with birth-time-marked (d — l)-dimensional maximal polytopes and 
makes the MNW-construction of Y(t, W) into a Markov process whose generator L is a 
clear modification of L as given in (5): 

LF(Y) = ( Yl ^ U K/}> S D " ^(^)]A(^) 

J\W] ^^IrZ^rrx 



/eCeiis(rnH) 



for F bounded and measurable on the space of birth time-marked tessellations of W. 



Consequently, writing p Y ( tA > w ^ — l ; . . . ,d — 1, for the birth-time-marked version of 



^Y{t\,w) ^ w ] iere g^jj A;-dimensional maximal polytope is marked with its birth time, by a 
straightforward modification of the proof of Theorem 2 we are led to 



Corollary 1 For all k = 1, . . . , d — 1 we have 

f Y^W) = {d _ k)2d -k-l j* I ^ P HT(sA,W) ^ ^ ^ 

6 Typical maximal polytope distributions 

We are now going to apply the results obtained in the last section to the stationary set-up, 
i.e. with A translation- invariant, to study the distribution of typical k- dimensional maximal 
polytopes of the STIT tessellation Y(tA) in M. d . Recall, that, in intuitive terms, the typical 
/c-dimensional maximal polytope of Y(tA) is what we get when we equiprobably choose 
a /c-dimensional maximal polytope of the tessellation out of a large observation window. 
This intuition can be made precise either by using Palm calculus or an ergodic approach, 
see [23, 24] for details. 
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Theorem 3 The distribution Q k of the typical k- dimensional maximal polytope ofY(tA) 
is given by 

/** rfs d_1 PHT(sA) , 

Jo ~l^^ k ' 

where Q^ HT<sA ' ) is the distribution of the typical k-face of the Poisson hyperplane tessellation 
with intensity measure sA. 

The preceding theorem can be rephrased by saying that the distribution of the typical 
A;-dimensional maximal polytope of a STIT tessellation is a mixture of suitable re-scalings 
of the distribution of the typical k- dimensional face of Poisson hyperplane tessellations, 
and that the mixing distribution is a beta-distribution on (0,t) with parameters d and 1. 



Proof of Theorem 3. To start, let (f k : MaxPolytopes fc — > R be a translation-invariant, 
non- negative measurable function for 1 < k < d — 1 and denote by Tp k {Y{tA)) the (possibly 
infinite) (/^-density of Y(tA) in the sense of [23, Chap. 4.1], i.e. 

lp k (Y(tA)) = lim 1 — E 

av 1 /eMaxPolytopes fc (y(*A,rWO) 

with W C R d some bounded convex set with positive and finite volume. The existence of 
this limit is guaranteed by Thm 4.1.3 ibidem. Using now Campbell's formula [23, Thm. 
3.1.2] and Theorem 2 from above, we obtain, possibly with both sides infinite, 



v ' /eMaxPolytopes fe (Y(tA,rWO) 



ds 



= (d-k)2 



d-k-l 



= (d-k)2 



d-k-l 



*1 

S 



lim — — - 

r^oo r d Vo\ d (W) 



Tp k (PRT(sA))ds, 



/eFaces fe (PHT(sA,rU')) 



(15) 



where PHT(sA) (or PHT(sA, rW)) is the stationary Poisson hyperplane tessellation with 
intensity measure sA (restricted to the window rW). We use the formula (15) with (p k = l 
and denote in this case Tp k (Y{tA)) =: \ k , referred to as the intensity of k-dimensional 
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maximal polytopes. Also denote by A^ HT< - sA ' 1 the intensity of fc-faces of the Poisson hyper- 
plane tessellation PHT(sA). By stationarity, there exists a constant G (0, oo), which is 
independent of s, such that A™ T( - sA ^ = CkS d . Thus, using (15), we obtain 



o 



A fc = (d - A;)2 d - fc - 1 / ic^d* = d —^2 d - k - 1 C k t d = d _±2 d - k -X^ tK) 

s ci (X 



(16) 



For a /c-dimensional polytope /, let m(f) be some translation-covariant selector (for exam- 
ple its Steiner point or the centre of gravity). Further, letting Q£ HT<sA ' ) be the distribution 
of the typical /c-face of PHT(sA) and be the distribution of the typical /c-dimensional 
maximal polytope of the STIT tessellation Y(tA), we obtain, again using (15), this time 
with tfk(f) := 1[/ — m W) G •], together with the definition of Qk and Q^ HT ^ A \ see [23, 
Eq. (4.8, 4.9)], 



td-k-l t 1 ^PHT^A^PHTfsA) 



A fe Q fe = (d - k)2 d ~ k ^ / -X^^Ql^'ds. (17) 



s 



Combining (17) with (16) we get, with Ck as above, 



t , , PHT(sA) r t , x PHT(sA) 



* ^^PHT( S A) ds = /•* ^ Q PHT( S A)^ 



which is the desired expression for Q k . □ 

It is interesting to note that formally marking the k- dimensional maximal polytopes with 
their birth-times and repeating the argument leading to Theorem 3 with Theorem 2 re- 
placed by its time-marked extension in Corollary 1 we obtain the birth time-marked ex- 
tension of Theorem 3: 

Corollary 2 The distribution Qk of the typical birth-time-marked k-dimensional maximal 
polytope ofY(tA) is given by 



f 

Jo 



ds d ~ 



t d 



ds. 



From these comparison theorems, mean values and even some distributional results for the 
typical dimensional maximal polytope of a stationary but anisotropic STIT tessellation 
can be concluded. To illustrate the general method, we exemplarily calculate at first 
the mean intrinsic volumes of the typical A;-dimensional maximal polytope of the STIT 
tessellation Y(tA). To neatly formulate them, let II be the associated zonoid of the Poisson 
hyperplane tessellation PHT(tA) in the sense of [23, Chap. 4.6]. 
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Corollary 3 The j-th intrinsic volume Vj (in the sense of integral geometry) of the typical 
k-dimensional maximal polytope Ik (with 0<j<k<d— 1) is given by 

W , M (ti) d wn) 

Proof of Corollary 3. Using Theorem 3, [23, Thm. 10.3.3] and the homogeneity of the 
intrinsic volumes, we get 

w.m - r dsd ~ i £^ (fn) .^ - v ^ r ^ 



o fflvou(fn) "~ © voi,(n)7 * (f) < 

(ti) d V d _ k (U) 

Q d-jvou(n)' 

which completes the proof. □ 

Note that in the isotropic case, i.e. when 1Z is the uniform distribution on Sd-i in the 
factorization (2), the zonoid II is a d- dimensional ball with radius proportional to t. More 
specifically, we have in this case 

mh) - * m^dirmvi 



(d-3)Kj\j) V r(f) y 

As a second example we turn now to the length distribution of the typical I-segment, which 
is nothing than the maximal polytope of dimension one. 

Corollary 4 The distribution of the length of the typical I-segment of a stationary and 
isotropic STIT tessellation with time parameter t > is a mixture of exponential distri- 
butions with parameter 71s. The mixing distribution is a beta- distribution on (0, t) with 
parameters d and 1. Its density with respect to the Lebesgue measure on (0, oo) is given by 



ft ds^ 1 d 

Pd(x) = J o llS e-^—ds = (7it)vm r(^+l,7ito), 



where T(-, •) is the lower incomplete Gamma- function and 71 = r(d/2)/(r(l/2)r((d + 
l)/2))- 

Proof of Corollary 4. This follows immediately from Theorem 3 and the well known 
fact that the length distribution of the typical edge of a stationary and isotropic Poisson 
hyperplane tessellation with intensity < s < t is an exponential distribution with param- 
eter 71 s. □ 
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In particular for d — 2 and d = 3 we have the densities 



p 2 (x) = ^(n 2 -(n 2 + 27rtx + 2t 2 x 2 )e-^, 

p 3 (x) = —r (48- (48 + 24:tx + 6t 2 x 2 + t 3 x 3 )e-^ tx ) . 

t 6 x V / 

The mean segment lengths are n/t in the planar case (see also [21]) and 3/t for d = 3 
(compare with [32]). Moreover, the variance of the length of the typical I-segment in the 
spatial case is given by 24/ 't 2 , which was not known before. In general, from the explicit 
length density formula it is easily seen that for the length of the typical I-segment only the 
moments of order 1 up to d — 1 exist. 
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